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Abstract 

We reconstruct the charmonium mesons J/tp , tpC^S) and Xc using a sample of 8.46 x 10^ BB 
events collected by the BABAR detector operating at e~^e~ center of mass energies near the T{4S) 
resonance. By measuring rates relative to the branching fraction of the J/tp , we obtain preliminary 
inclusive B branching fractions of (0.25 ± 0.02 ± 0.02)% to the ^'(25') and (0.39 ± 0.04 ± 0.04)% to 
the Xci, and set a 90% confidence level limit of 0.24% on decays through the Xc2- 
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1 Introduction 



Observation of the charmonium mesons J/ij: , ip{2S) and Xc is a critical component of the mea- 
surement of CP violation in B decays at BABAR and other facilities. Measurements of the in- 
clusive branching fractions of B mesons to these states are made using the decays J/ip — > , 
■0(25) — > ■0(25) -K'^TT'J/ijj and Xc — *■ 7-^/0 1 with J/0 — > i'^l" in the latter two cases; 

the leptons £ may be either electrons or muons. The production rates of the ip{2S) and Xc mesons 
are measured relative to the J/ij: in order to eliminate common systematic errors in tracking and 
particle identification. 

2 The BaBAR detector and dataset 

The data used in this analysis were collected with the BABMi detector at the PEP-II storage ring, 
which collides 9.0 GeV electrons with 3.1 GeV positrons. A luminosity of 7.7 fh~^ was collected on 
the T(4S') resonance, y/s = 10.58 GeV ( "on-resonance" ) , with an additional 1.2 fh~^ collected at 
y/s = 10.54 GeV, below the threshold for BB production ("off-resonance"). The off-resonance data 
set is used to statistically subtract the non-BB component of signals. 

The BABAR detector is described elsewhere [^. The tracking system, which is used for pattern 
recognition and reconstruction of the momenta of charged particles, consists of two sub-detectors. 
The inner is a five-layer double sided silicon vertex tracker (SVT) which gives precision spatial 
information for all charged particles and is the primary detection device for low momentum charged 
particles. It is surrounded by a 40-layer drift chamber (DCH), which provides measurements of 
track momenta. 

The primary sources of information used in lepton identification are the calorimeter and the 
muon system. The DCH also provides dE/dx information which is used in the identification 
of electrons. The calorimeter is constructed from Csl crystals and, in addition to distinguishing 
electrons, hadrons and muons on the basis of their energy depositions, is used to locate and measure 
photons. The muon system, known as the Instrumented Flux Return (IFR) in BABAR, consists of 
resistive plate chambers interleaved with iron plates. There are 19 layers in the central region of 
the detector, 18 in the endcap regions. 

3 Event selection 

The event selection criteria are designed to select BB decays with high efficiency while suppressing 
lepton pairs, two-photon events, and interactions with the beam pipe or residual gas in the beam 
line. We require the event to have: 

• at least four charged tracks in the active volume of the tracking system; 

• total energy (charged tracks plus photons) greater than 5 GeV; 

• a primary vertex within 5 mm of the beam interaction region; and 

• a normalized second-order Fox- Wolfram moment [|[ less than 0.7. 

A total of (8.46 lb 0.14) x 10^ BB events satisfy these criteria. This value is obtained by scaling the 
number of events selected in the off resonance sample by the ratio of muon pairs (luminosities) and 
subtracting it from the total number of events selected in the on-resonance sample. The uncertainty 
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Figure 1: Mass distribution of J/i/j candidates to (a) electron pairs and (b) muon pairs after the 
continuum subtraction. 



is dominated by variations in the ratio of selected events to muon pairs in the off-resonance running 
periods throughout the data collection period. 



4 Reconstruction of charmonium mesons 
4.1 Reconstruction of J/tp 

J/ip candidates are formed from pairs of charged tracks in the active fiducial volume of the tracking 
system and calorimeter. Tracks are required to pass within 1.5 cm of the beam line and to include 
information from the DCH. 

Both tracks are required to satify lepton-identification criteria. Electrons must have an energy 
deposition in the calorimeter of at least 75% of that expected from the track's momentum and 
a shape consistent with an electromagnetic shower. Muons must have an energy deposition in 
the calorimeter of less than 0.5 GeV — approximately twice the expectation for a minimum ionizing 
track — and traverse an amount of material that is within two interaction lengths of that expected 
for a muon of the measured momentum. The pattern of energy deposition in the IFR must be 
narrow, perpendicular to the flight path of the track, and the location of the hits must have an 
acceptable when fit to the trajectory extrapolated from the DCH. 

The two tracks are fit to a single vertex, if possible, before calculating the mass of the candidate 
J/ip . Otherwise, the mass is derived from a simple sum of their four-vectors. 

Histograms of mass distributions are accumulated for both the on- and off-resonance data sets. 
These are then subtracted, after weighting for the difference in luminosities, to give "continuum 
subtracted" distributions. Figure |l] shows the resulting mass distributions for J/-0 — > e~^e~ and 
J/ip fJ-^fJ'^ candidates. 

The number of J/'0 mesons is extracted from the mass distribution using a probability distribu- 
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tion function (p.d.f.) derived from a simulation that includes final state radiation and bremsstrahlung, 
The mass resolution of the detector simulation is increased from 11.0 to 13.3MeV/c^ to match the 
observed value. 

A systematic error on the fit to the e~^e~ final state is obtained by varying by ±10% the number 
of J/ip mesons in which an electron undergoes bremsstrahlung in the simulation. This variation, 
which reflects possible differences between the detector simulation and reality, translates to a 3.5% 
systematic error. 

Systematic errors for both final states are obtained by modifying the selection criteria in two 
ways. First, the momentum of the J/ip candidate in the T{4S) center of mass (CM) system, p*, is 
required to be less than 2.0GeV/c; all J/V' mesons from B decays should pass. Second, the number 
of charged tracks required is increased to five in order to reject remaining lepton pair events. These 
studies produce systematic errors of 1.3% for e~^e~ and 1.6% for fi'^n". 

We find 4920 ± 100 ± 180 J/ip e+e" and 5490 ± 90 ± 90 J/V' ^ /i+zx" events with mass 
greater than 2.6GeV/c^. 

4.2 Reconstruction of 'ip(2S) 

The reconstruction of the ip{2S) through its decay to e'^e'^ or is similar to that of the J/-0 . In 

order to reduce the relatively higher background rates, additional selection criteria are applied to 
the ip{2S) candidate. First, its momentum in the CM system (p*) must be less than 1.6 GeV/c, true 
for all ip{2S) mesons from B decays. Second, the two tracks must be consistent with coming from 
a common vertex with a probability greater than 1%. Finally, more stringent particle identification 
criteria are applied. Muon candidates must penetrate to within one interaction length of the 
expected amount of material and at least one electron candidate must have an energy deposition 
between 0.88 and 1.30 times its momentum. 

The number of signal events is extracted from the mass distribution using a p.d.f. from simula- 
tion, as for the J/ip . We do not perform a continuum subtraction in this case, as the requirement 
on p* is expected to remove most backgrounds from the continuum. Variation of the amount 
of bremsstrahlung in the detector simulation leads to a 1.7% systematic error on the number of 
ip{2S) — > e^e~ events, although this error is negligible compared to the statistical errors in this 
case. We find 131 ± 29 ± 2 ^j{2S) e+e~ and 125 ± 19 ij{2S) fi+n' events (Fig. |). 

Due to their lower momentum, the pion candidates used to reconstruct ip{2S) — > 7r"'"7r~ J/-0 are 
not required to have DCH information. Pairs of pions with mass between 0.45 and 0.59 GeV/c^ 
are combined with J/tp candidates with mass between 2.60 and 3.13GeV/c2 for e+e" and 3.06 
and 3.13GeV/c^ (—2.8 to +2.5 standard deviations) for /i^/x". The wider window for e^e^ is 
necessitated by the larger radiative tail in this final state. The above requirements on p* and on 
vertex quality are applied in this final state as well. 

To reduce the impact of the radiative tail and the mass resolution of the J/tp candidate, the 
number of ip(2S) events is extracted from a fit to the distribution of the mass difference between 
the il)(2S) and the J/ip (Fig. |3|), with a Gaussian distribution for the signal and a Chebychev 
polynomial for background. In candidates with J/ip — > e^e~ , we obtain 126 it 44 ip{2S) events, 
while for J/ip fj,'^fj,~, we obtain 162 it 23. 

4.3 Reconstruction of Xc 

Xc mesons are reconstructed through the radiative decay Xc iJ/i' ■ The selected photon candi- 
dates must be reconstructed within the same fiducial region used for leptons and have an energy 
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Figure 2: Mass distribution of tp{2S) candidates to (a) electron pairs and (b) muon pairs. 




0.5 



0.55 



J/v|/ e^e' 

_l I I I I I I I I I l_ 



^ I I I I I 



b) 



^ I I I I I 



B\RiR 




_l I I I I I I L_ 



J/V|/ i^V" 

J I I I I I I I I l_ 



0.6 



0.65 



0.7 



0.5 



0.55 



0.6 



0.65 



0.7 



M(e+e-7U+7U-)-M(e+e-) (GeV/c ^ 



M(^l^^l"7u+7u-)-Mal^^l") (Oev/c ) 



Figure 3: Mass difference between the 'tp{2S) and J/ip candidates, where tp{2S) — >■ tt+tt J/tp and 
(a) J/ip — > e"'"e~ and (b) J/i/j IJ'^IJ'~- 
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between 0.20 and 0.65 GeV in the CM system. The shape of the energy deposition must be consis- 
tent with an electromagnetic shower. To reject photons from nearby hadronic showers, the photon 
must be at least 9° from the nearest charged track in the calorimeter. The substantial combina- 
torial background from vr'^ decays is reduced by rejecting any photon that, when combined with 
any other photon, forms a tt^ candidate with mass between 0.109 and 0.147 GeV/ c'^ (—4.0 to +2.0 
standard deviations about the vr'^ mass). 

In selecting the J/'0, muon candidates must satisfy the more stringent requirement that they 
penetrate to within one interaction length of the expected amount of material. J/^l^ candidates with 
mass between 3.05 and 3.14 GeV/ are combined with photon candidates to form Xc candidates. 

As for the ip{2S)^ the number of events is extracted by fitting the distribution of the mass 
difference between the Xc and J/ijj candidates (Fig. The resolution is dominated by the angular 
and energy resolutions for the photon, and includes a low-side tail due to preshowering before the 
calorimeter. Xc signals are extracted from a fit with a "Crystal Ball" function, which includes this 
feature. The parameter that describes the tail is fixed to the value determined from a complete 
detector simulation |^]. 

A 4.2% systematic error due to a possible mismatch between the fitting function and the true 
distribution is estimated from fits to simulated data. An additional systematic error is evaluated 
by requiring p* < 2.0GeV/c for the J/-0 candidate. The resulting variation in the signal, 9.3% for 
e~^e~ and 2.0% for fj,^^~, although not inconsistent with a statistical fluctuation, is taken as a 
systematic error. 

We simultaneously fit for a Xci and a possible Xc2 component. The only additional free parame- 
ter is the number of Xc2 candidates: the mass difference between the Xc2 and the Xci is fixed to the 
Particle Data Group value and the resolution and tail parameters are fixed to those of the XciQ- 

We find 129 it 26 it 13 Xci and 3 it 21 Xc2 events in candidates in which J/ip — > e~^e~ and 
204 lb 47 lb 12 Xci and 47 ib 21 Xc2 in candidates in which J/tl^ ^u"*"//". 

5 Calculation of branching fractions 
5.1 B^J/ipX 

The number of J/tp mesons produced in our data sample is related to the number extracted from the 
fit by a correction for the reconstruction efficiency and by the branching fraction for J/il^ — > 
We average the published J/ip e~^e~ and J/ip /^^ branching fractions, assuming a common 
systematic error, to get (5.91 ib 0.10)%[|]. 

A complete detector simulation is used to extract the probability that both leptons are recon- 
structed as tracks of sufficient quality within the fiducial volume, with a combined mass greater 
than 2.6 GeV/ (the region used in the fit to the data). 

The accuracy of the detector simulation in reproducing tracking efficiency has been studied in 

^ K- vr^TT+vr decays and in r pairs in which one r decays to a single charged track and the 
other decays to three charged track From these studies, a uncorrelated systematic error of 2.5% 
per track is assigned to the overall efficiency. 

The probability of the leptons being identified correctly is determined using "control samples" — 
sets of tracks in which the nature of the particle is known without the use of the particle iden- 
tification systems. For both electrons and muons, we use lepton pairs and radiative lepton pairs 
produced in two-photon events. A track is included in the control sample if the kinematics of the 
event are consistent with expectations and if the other track in the event passes strict particle 
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Figure 4: Mass difference between tfie Xc and J/if) candidates, where Xc iJ/i^ a^nd (a) J/ip 
e+e~ and (b) J/ijj — ■ 

identification requirements. Tfic control sample tracks are then used to evaluate the efficiency of 
the selection criteria as a function of momentum and location in the detector. A value for the 
branching fraction B J/ip X will be determined once this study has been completed. 

5.2 B ijj(2S)X 

To minimize systematic errors from tracking and lepton identification, we calculate the ratio of 

'0(25') and J/^p inclusive branching fractions using the ratio of observed ^{2S) to J/ip mesons. 
For the purposes of this calculation, we refit the J/tp signal after applying the more stringent 
particle identification criteria used in reconstructing the ip{2S). Simulation studies indicate that 
the probability for the il^{2S) — £'^£~ leptons to be reconstructed and pass the identification criteria 
is equal to that for the J/ip £'^i~ decay. 

The ratio of iIj{2S) to J/ip inclusive branching fractions is equal to the ratio of observed mesons, 
corrected for acceptance, the branching fractions J/-0 t^l~ and ip{2S) and a correction 

for the efficiency of the vertexing requirement. The vertexing efficiency is found by simulation to 
be 0.941 with a 2% systematic error, which is determined from vertexing studies of the J/tp . 

For the branching fraction of 1^(23) i'^i^ we average the two Particle Data Group values 
and obtain (0.90 zt 0.12)%. This is the dominant systematic error. 

The ratio of ip{2S) to J/ijj inclusive branching fractions is calculated from the tt+tt" J/^ final 
state data in a similar fashion. It is equal to the number of reconstructed ip(2S) mesons divided 
by the number of J/ip mesons reconstructed within the tighter mass windows used in this analysis, 
with corrections for the pion reconstruction probability, the vertexing efficiency and the ijj{2S) 

J/if) branching fraction. The probability for the two pions to be reconstructed within the 
fiducial requirements and the mass window is 0.524. The efficiency of the requirement for a good 
vertex is found by simulation to be 0.89 with a 4% systematic error. The uncertainty in the 
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published branching fraction for ip{2S) tt^tt^ J/ip , (31.0 it 2.8)%, is the largest component of the 
systematic error. 

The two decay modes give consistent values for the ratio of inclusive branching fractions B{B ^ 
'ip{2S)X)/B{B J/ipX) and are combined to give an overall value in Table |l|. The ratio is 
multiplied by the Particle Data Group value for B{B — > J/-0X), (1.15 it 0.06)%, to obtain the 
inclusive iIj{2S) branching fraction. 

5.3 B XcX 

The ratio of Xci to J/ip inclusive branching fractions is computed from the number of reconstructed 
ip{2S) mesons divided by the number of J/ip mesons reconstructed within the tighter mass windows 
and tighter particle identification criteria used in this analysis, the published Xci — > lJ/''P branching 
fraction, and a correction for the photon reconstruction probability. 

The efficiency for the photon reconstruction is calculated from the simulation to be 0.550. A 
systematic error of 1.6% is obtained by varying the calorimeter simulation resolution and energy 
scales over bounds determined by studies of vr'' reconstruction. The branching fraction for Xci 
^J/ip, (27.3 lb 1.6)%, contributes a 5.8% systematic error. The systematic errors on the fit have 
been described earlier. 

To obtain a limit on the production of the Xc2, the fit results from Fig. ^ are recast in terms 
of the ratio of Xc2 to J/V' inclusive branching fractions. The likelihood function for each fit is then 
normalized to unity for N^.^^ > and multiplied together to give the combined likelihood. The 
90% confidence level is taken to be the ratio of branching fractions that bounds 90% of the area of 
the likelihood function. 

The statistical errors from the fit lead to a 90% confidence limit on B{B — > Xc2)/B{B J/ip) of 
0.186. However, there is an 8% systematic error in converting the fits into limits on Xc2 production 
due to the uncertainty in the published branching fraction for Xc2 iJ/i^- We include this 
uncertainty by increasing the limit by 1.28 times this error, giving a limit of 0.205. 

The results for Xci and Xc2 are summarized and presented as inclusive B branching fractions 
in Table 0. 

Table 1: Summary of inclusive B branching ratios measured with respect to i3 — > J/ipX. The total 
systematic error in the e'^e~ and final states is the sum in quadrature of the final-state- 

specific values and the systematic error common to both. The results are combined and multiplied 
by the PDG value for B{B — > J/ipX) to obtain inclusive branching fractions (%). 





Branching ratio relative to J/V' 


Branching fraction (%) 




e+e 

Stat Sys 


Stat Sys 


Common 

Sys (%) 


Combined 
Stat Sys 


Xci jJ/i^ 


0.28 0.06 0.03 


0.38 0.05 0.02 


7.3 


0.39 0.04 0.04 


Xc2 90% CL 


0.16 


0.28 


8.1 


0.24 


ij{2S) £+£- 
^{2S) ^+7r- J/V- 
Combined 'ip{2S) 


0.22 0.05 0.01 
0.18 0.06 0.01 


0.23 0.04 0.01 
0.22 0.03 0.00 


13.6 
9.9 


0.26 0.03 0.04 
0.24 0.03 0.03 
0.25 0.02 0.02 
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6 Summary 



The charmonium states J/'tp , i/'(2S') and Xc have been reconstructed in B meson decays. Prehminary 
measurements of the ratio of branching fractions B{B tlj{2S)X)/B{B J/ipX) and B{B 
XciX)/B{B J/jpX) have been presented and a hmit set on the ratio for Xc2- Using the Particle 
Data Group value for the inclusive branching fraction for J/ijj, these values have been converted 
into inclusive B branching fractions. 

A direct value for the J/ip branching fraction will be available upon completion of particle 
identification studies. Several of these results are statistically limited; more precise values will 
be available based on data that are currently being recorded. The new data will also permit 
measurements of the p* and helicity distributions for these charmonium mesons. 
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